Baryon and meson screening masses 
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In a strongly-coupled quark-gluon plasma, collective excitations of gluons and quarks should dom- 
inate over the excitation of individual quasi-free gluon and quark modes. To explore this possibility, 
we computed screening masses for ground-state light-quark mesons and baryons at leading-order in a 
symmetry-preserving truncation scheme for the Dyson-Schwinger equations using a confining formu- 
lation of a contact-interaction at nonzero temperature. Meson screening masses are obtained from 
Bethe-Salpeter equations; and baryon analogues from a novel construction of the Faddeev equation, 
which employs an improved quark-exchange approximation in the kernel. Our treatment implements 
a deconfinement transition that is coincident with chiral symmetry restoration in the chiral limit, 
when both transitions are second order. Despite deconfinement, in all T = bound-state channels, 
strong correlations persist above the critical temperature, T > Tc, and, in the spectrum defined 
by the associated screening masses, degeneracy between parity-partner correlations is apparent for 
T > 1.3 Tc. Notwithstanding these results, there are reasons (including Golberger-Treiman relations) 
to suppose that the inertial masses of light-quark bound-states, when they may be defined, vanish at 
the deconfinement temperature; and that this is a signal of bound-state dissolution. Where a sensi- 
ble comparison is possible, our predictions are consistent with results from contemporary numerical 
simulations of lattice-regularised QCD. 

PACS numbers: 12.38.Aw, ll.30.Rd, 12.38.Lg, 24.85. +p 



I. INTRODUCTION 

It is generally held [l[ that a strongly-coupled quark- 
gluon plasma (sQGP) has been produced at the relativis- 
tic heavy ion collider (RHIC) and that this substance 
behaves as a nearly perfect fluid on a domain of temper- 
ature that extends above that required for its creation: 
If such is the case, then collective 
excitations of gluons and quarks should dominate within 
this domain over the excitation of individual quasi-free 
gluon and quark modes. It has long been knovifn that 
one may explore this possibility by studying the screen- 
ing masses of hadrons above Tc [3| , since the long-range 
structure of a plasma determines quantities such as its 
equation of state and transport properties. It is note- 
worthy, however, that sufficiently removed from T = 0, 
the screening masses in a given channel do not have a 
simple connection with the inertial masses of what were 
the channel's bound-states at T = 0. 

In the calculation of masses (screening or inertial), as 
with most other applications, in order to arrive at robust 
conclusions, one should employ an approach that pre- 
serves the symmetries of QCD, is simultaneously appli- 
cable to both mesons and baryons, and has already been 
applied with success to a diverse array of observables. 
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The numerical simulation of lattice-regularised QCD is 
one such approach, from which some results, relevant 
herein, are described, e.g., in Refs. It is found that 

correlations persist in hadron channels above Tc and that 
parity-partner screening masses become equal. 

The temperature dependence of hadron masses may 
also be explored in the continuum using models with 
an arguable foundation in QCD. Unfortunately, few such 
models are simultaneously applicable to studies of mesons 
and baryons; and amongst those that arc, different mod- 
els and treatments lead to different conclusions. 

For example, in a formulation of the Global Colour- 
symmetry Model that treats mesons and baryons dif- 
ferently, it is found Q that the nucleon's mass vanishes 
and the associated radius diverges at the chiral symme- 
try restoration temperature (Tc). Alternatively, a scalar- 
diquark-only Faddeev equation treatment of the nucleon, 
within a Nambu-Jona-Lasinio model framework Q that 
docs not express confinement, indicates [l^ that whilst 
the diquark correlation can survive above the tempera- 
ture for chiral symmetry restoration, the nucleon docs 
not: its mass falls with increasing temperature and, be- 
low Tc, it lies above the quark-|-scalar-diquark breakup 
threshold. One might nevertheless argue that there is 
some qualitative similarity between these outcomes. 

On the other hand, an analysis of the nucleon using 
thermal finite-energy sum rules produces [ll| a nucleon 
mass that rises rapidly in the neighbourhood of Tc, owing, 
it is suggested, to a vanishing probability for three quarks 
to form a correlation above this temperature. In this 
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framework, results for mesons are similar [12|. 

These results are all interpreted as statements about 
inertial masses; and deconfinement is claimed in each 
framework, although the signal is markedly different. 
Plainly, therefore, results obtained in the computation 
of inertial masses at ttT > Aqcd arc very sensitive to 
assumptions; and such discordant outcomes are unsat- 
isfactory. This prompts us to re-examine the evolution 
of hadron masses with temperature in the continuum. 
We will, however, focus primarily on the computation of 
screening masses. These static quantities are only indi- 
rectly sensitive to the real-time scattering processes that 
may complicate the evaluation of inertial masses, which, 
in-medium, are all likely to become complex numbers: 
M m{T)-iuj{T) [H]. Finally, however, and m princi- 
ple, upon computation of all screening masses in a given 
channel, the inertial masses may be reconstructed. 

We choose to employ QCD's Dyson- Schwingcr equa- 
tions (DSEs) in our analysis. The DSEs have long been 
used to explore the response of hadron phenomena to 
nonzero temperature and density (l^ [isj . They are an 
excellent tool for use in our study because they provide 
a Poincare-covariant framework that p^| - [l8j : is capa- 
ble of simultaneously implementing light-quark confine- 
ment and expressing dynamical chiral symmetry break- 
ing (DCSB); admits a symmetry-preserving truncation 
scheme; and provides a unified and uniform treatment of 
mesons and baryons 

Here it is worth commenting further on the manner by 
which the study of mesons and baryons is unified. Both 
are treated as continuum bound-state problems with ker- 
nels built from dressed-gluon and dressed-quark propa- 
gators that express the dynamical generation of mass: 
mesons via Bethe-Salpeter equations (BSEs) and baryons 
via Faddeev equations. In this way, they are completely 
equivalent. Formulation of the Faddeev equations is sim- 
plified by capitalising on the established importance of di- 
quark correlations [l9l.[20|. This is not unusual. However, 
in extending the Faddeev equations to nonzero temper- 
ature herein, we retain axial-vector diquark correlations. 
They are crucial at zero temperature [HI, [13]; e.g., they 
provide significant attraction in the nucleon bound-state 
and the A-resonance is inaccessible without them. There 
is no reason to expect axial- vector diquarks to be less im- 
portant at r 7^ 0; and their omission in earlier studies 
undermines the reliability of those analyses. 

In employing the DSE approach herein, we will exploit 
a symmetry-preserving treatment of a vector x vector con- 
tact interaction because it produces T = results for 
hadron static properties that are not realistically distin- 
guishable from those obtained by using more sophisti- 
cated kernels. This is demonstrated in Refs. [2ll [23l - |28| . 

In Secin] we explain our confining, symmetry- 
preserving treatment of the contact interaction and its 
extension to nonzero temperature. This Section covers 
the gap, Bethe-Salpeter and Faddeev equations. It is 
augmented by two appendices, which detail the deriva- 
tion of the Faddeev equations. We highlight at the 



outset that our Faddeev equations do not express the 
full complexity that arises in-medium. Notwithstanding 
this, they are practical simplifications which should at 
least yield qualitatively and semi-quantitatively reliable 
insight into the behaviour of nucleon and A screening 
masses. 

Our formulation of the contact interaction implements 
confinement, following the notions of Ref. [2^. This dis- 
tinguishes it from most previous simultaneous, covariant, 
continuum treatments of mesons and baryons in-medium. 
On the other hand, we expect QCD to exhibit deconfine- 
ment at some T = > 0. Consequently, a dynamical 
mechanism should be incorporated the allows for decon- 
finement. In a related context, this was considered in 
Ref. [lOl- Our different approach is described in Sec. lIIII 

Our results are described in Sec. lIVI They range 
from an analysis and illustration of chiral symmetry 
restoration and deconfinement to a description of the 
T-dependence of screening masses of ground-state light- 
quark mesons, diquark correlations and baryons. A key 
result is that correlations typically persist in hadron 
channels above the temperature associated with chiral 
symmetry restoration and deconfinement. As we explain, 
however, this is not a statement that true bound-states 
persist. 

In Sec.|Vl following Ref. [IH, we explore the implica- 
tions of chiral-symmetry restoration and deconfinement 
for quark- and nucleon-level Goldberger-Treiman rela- 
tions; and therefrom, the nucleon's inertial mass. Sec- 
tion |VI] presents a summary and perspective. 



II. CONTACT INTERACTION AT NONZERO 
TEMPERATURE 

The formulation of DSEs at nonzero temperature is 
described in Refs. so here we proceed directly with 

a specific discussion of the contact interaction described 
in App.S 



A. Dressed quark propagator 

The T 7^ dressed-quark propagator is obtained from 
the following gap equation: 

Jl,dq 

where m = m„ = is the light-quark current-mass, 
= Tj:Z.^Jd^q/{2nf and = (2n + 1)^T is 
the fermion Matsubara frequency. Equation ([T]) is the 
rainbow approximation to the gap equation, which is 
the leading term in the systematic, symmetry-preserving 
truncation scheme of Ref. [H, . 
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The solution of Eq. ([T]) is momentum-independent, as 
it was for T = 0; viz., 



S ^ (p, iu„) = 17 • p + i^iiOn + M , 

where the dressed-quark mass is determined from 

, 167raiR f Ud 

Ad = ni + ■ 



(2) 



(3) 



with si = + Lof . The integral is linearly divergent and 
may be regularised using the T = procedure; namely, 
we write 12911 



1 



Si 



AP 



Si 



M2 



(4) 



(5) 



where Ti^^uv are, respectively, infrared and ultraviolet reg- 
ulators. Since the interaction in Eq. ()A1|) does not define 
a renormalisable theory, then Auv '■— I/tuv cannot be 
removed but instead plays a dynamical role, setting the 
scale of all dimensioned quantities. Using Eq. ([5]), Eq. ([3]) 
becomes 



M = 



M 



4aiR. 



(6) 



where C™{Ad'^;T) is defined in Eq. ((X4| and the param- 
eters in Table H 

It is apparent from the rightmost expression in Eq. ([5]) 
that a finite value of Th- = 1 / Air implements confinement 
by ensuring the absence of quark production thresholds 
in all processes [13, [130 This is appropriate for studies of 
T = phenomena. However, we expect QCD to exhibit 
deconfinement at some T = T^; > 0, whereat the produc- 
tion thresholds reappear, as illustrated in Refs. [4^ |4^ . 
Therefore, we subsequently introduce a dynamical mech- 
anism that makes rir temperature-dependent. 

In the gap equation above, and in the bound-state 
equations to follow, we omit the temperature-induced 
separation of the gluon propagator dressing into trans- 
verse and longitudinal parts. Instead, we assume the 
interaction dressing is frozen at its T = form. This is 



a defect that our study shares with all other continuum 
analyses of bound-state screening masses. In our judge- 
ment, however, it is not a crippling weakness. Tempera- 
ture does affect the nature of gluon dressing but, on the 
domain of concern to us; i.e., T < 2Tc, these effects are 
modest 



B. Mesons 

1. Screening masses 

A strength of the DSE framework is its ability to treat 
mesons and baryons on an equal footing, both at zero 
and nonzero temperature. We illustrate and exploit that 
capacity herein. 

At leading order in the symmetry-preserving trunca- 
tion scheme of Ref. [11] , one considers the homogeneous 
rainbow-ladder BSE for mesons; namely. 



l.dt 



(7) 

where S is obtained from Eq. ([T]) and Qo = {Q, 0} is the 
total momentum entering the amplitude. (For our imme- 
diate purposes, it is necessary to focus only on the me- 
son's zeroth Matsubara frequency.) The rainbow-ladder 
truncation is known to provide reliable results for the 
T = properties of vector and flavour nonsinglet pseu- 
doscalar mesons [13, 113 • Equation ([7]) is an eigenvalue 
problem: it has a solution for Qg = —rnjj, where nin is 
the mass of any one of the bound-state's in this channel 
at T = 0, owing to 0(4)-invariance, and the associated 
screening mass for T > 0. 

We reiterate here that, sufficiently removed from T = 
0; namely, for ttT > M(T = 0), a particular screening 
mass in a given channel has no simple connection with the 
inertial mass of any of the channel's bound-states at T = 
0. Indeed, a given channel supports many bound-states 
at T = and hence there will necessarily be at least as 
many screening masses. Moreover, traceable to each T = 
bound-state, there is potentially a different screening 
mass for each one of the enumerable infinity of (meson 
or baryon) Matsubara modes. In order to reconstruct 
a real-time Green function, from which an inertial mass 
may be determined, one must: calculate each such mass; 
compute a Fourier transform of the form 



The potential between infinitely-heavy quarks measured in nu- 
merical simulations of quenched lattice-regularised QCD - the 
so-called static potential - is not related in any known manner 
to the question of confinement in the real world, in which light 
quarks are ubiquitous. It is a basic feature of QCD that light- 
particle creation and annihilation effects are essentially nonpcr- 
turbative and therefore it is impossible in principle to compute 
a potential between two light quarks l34l ISal . Confinement may, 
instead, be related to the analytic properties of QCD's propaga- 
tors and vertices |36|442(I . 



oo 

S{t) = tY^ e— "-5(i.„) , (8) 

n— — oo 

where {vn} are boson or fermion Matsubara frequencies 
and 5 is a thermal Schwinger function; complete an an- 
alytic continuation via r — > r + if and the limit t — > 0"*" ; 
and finally arrange appropriate step-function-weighted 
combinations of the result. 
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2. Pseudoscalar- and vector-mesons 

It was shown in Ref. [1^ that a pseudoscalar meson 
must possess components in its Bethe-Salpeter ampH- 
tude that may be described as pseudovector in charac- 
ter. These components play a critical role in the T — 
physics of pseudoscalar mesons 0, [13, H^, Hi], . This 
means herein that Eq. ([7]) supports a solution in the pseu- 
doscalar channel of the form: 

r.(Qo) = il^EAQo) + ^757 • QoFAQo) ■ (9) 

In the vector channel there are two distinct components 
at T ^ 0: 



rp(Qo) 



(10) 



where 7j_ = PijiQohj = (^-y - QiQj/\Q\'^h], = 
1,2,3. In a symmetry preserving treatment of the in- 
teraction in Eq. (|A1[) , there can be no dependence on a 
relative momentum in either case. 

Explicit forms for the BSEs of the tt- and p-mesons 
are readily obtained following the procedures described 
in Refs. [25|,[23|. That for the pion is an obvious analogue 
of Eqs. (31)-(35) in Ref. [H. Herein, as r"(cr) is replaced 
by r^(^;r), so 



C, M^Ci (qT) 



r"(qT). 



(11) 



The argument of these functions is given by (d = 1 — a) 



(12) 



Note that we follow the interaction definition of Ref. [27[ , 
so that I/tjiq in Ref. [Isj becomes iTTam/rriQ herein. 
(See Eq. (20) in Ref. [131 ■) For brevity, we will sometimes 
write 



1 



Ana 



IR 



(13) 



Given these observations, one can readily express the 
pseudoscalar BSE: 



Et, (Qo ) 


4aiR 






' E4Qo)' 










_ F^Qa) _ 



where 



^E 

Aff 



C da [C^'i^iA'P.a.-mD-T) 
Jo 

+2aamlCi{'^iAd'^,a,-ml);T)] , (15a) 
daC'"{<;{Ad^,a,-ml);T), (15b) 



(15d) 



hl^ J dad^{^{M\a,-ml)-T), (15c) 
-IKfe- 



In deriving Eqs. (|T4)) . (|T5]) . one must use the following 
identity: 



0= da 
Jo 



c'"(qr) + a"(qr) + 5(.'"(qr) 



(16) 



where 5(.'"(<?;T) is defined in Eq.ESJ. As the T ^ 
generalisation of Eq. (20) in Ref. [25[, Eq. is neces- 
sary and sufficient to guarantee the vector and axial- 
vector Ward-Green- Takahashi identities [50l - [52j . N.B. 
^'"(cr;T 0) = 0, as shown in connection with 

Eqs.dMl), ([Sl- 
it is necessary to employ the canonically normalised 
Bethe-Salpeter amplitude in the computation of observ- 
ables. For the pion, that amplitude satisfies 



_d_ 

dQl' 



1= -TT^^AK^Qo) 



(17) 



K=Qo 



where 

n.(A', go) 



-etro / T^{~K)S{q + Q,cJi)r^{K)S{q,iJi). (18) 

Jl,dq 

The BSE for the p-meson's longitudinal component is 



(19) 



where 



4q;ir 



da 



aazCi {'^{A'l^ ,a, z);T) 



+5(.'"(?(M2,a,z);r) 



(20) 



The canonical normalisation condition for the longitudi- 
nal amplitude is 



(21) 



(14) Z 



The BSE for the perpendicular component of the p- 
meson is almost identical to Eq. ([TS]) except that one 
omits ^™ in mapping %^ ■ The equation thus 

obtained is just Eq. (36) in Ref. [2^ with the replacement 

C\ {s) ~^ Ci {s\T\ The canonical normalisation condi- 
tion is simply Eq. (PT|) with ||— >-L in the obvious places. 
At zero temperature the equations for the longitudinal 
and transverse components are naturally identical, an 
outcome that is realised because %^^(<,;T — >■ 0) =0. 



3. Scalar and pseudovector channels 

The large splitting between parity partners is a striking 
feature ofQCD's spectrum. It is discussed at length in 
Refs. [5l,[55, which showed that nonperturbative DCSB- 
corrections to the rainbow-ladder truncation generate 
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a large spin-orbit repulsion and are responsible for the 
splitting. Informed by those analyses, Refs. H^, mod- 
ified the rainbow-ladder BSEs in the scalar and pseu- 
dovector channels, including a single, common coupling 
parameter g^Q whose presence simulates the repulsive ef- 
fect. The value 



T=0 

9so 



0.24 



(22) 



reproduces the experimental value for the ai-p splitting. 
It is noteworthy that the shift in is accompanied 
by an increase of rria, the new value of which matches 
an estimate for the ^(/-component (dressed-quark core) 
of the cr-meson obtained using unitarised chiral pertur- 
bation theory [55, , 5&| . 

We emulate Refs. m,!!^ by including gso in our BSEs 
for the scalar and pseudovector channels. However, in 
anticipation that chiral symmetry is restored above some 
critical temperature, we enable the parameter's strength 
to track that of DCSB; viz., we use 



5io->5io(r) = l 



M{T) 



(1 



(23) 



where M{T) is the T-dependent dressed-quark mass ob- 
tained from Eq. ([3]) . 

With a symmetry-preserving regularisation of the con- 
tact interaction, the scalar meson Bethe-Salpeter ampli- 
tude takes the simple form 



r,(Qo) = lD£^a(Qo) 
The screening mass is obtained from 
4q!ir 



1 = 



(24) 



(25) 



where 

%_^{z)=gl^{T) \\a\C^(s{M\a,z)-T) 
Ja 

-2C(<^(Af^a,z);^)]; 
and the amplitude is canonically normalised via 
1 3 d 



El 27r2 d;. 



(26) 



(27) 



Akin to the p-meson, at nonzero temperature the 
Bcthe-Salpeter amplitude for the ai channel has the form 



7574-B^i(Qo) 

757±i?i(Qo) 



(28) 



The screening mass of the transverse component is ob- 
tained from 



= 1 



(29) 



where 



gloiT) 



4aiR 



daCr{^{M'^,a.z)-T). 



(30) 



The canonical normalisation condition is 



(31) 



The BSE for the longitudinal component is modified 
similarly to that of the p, Eq. ([20]) : 



= 1 



(32) 



where 



\,{z)^glo{T) 



4aiR 
37rm2, JO 



da [cr(<?(M2,a,z);T) 



+!^-(?(M2,a,z);r)] ; (33) 
and the canonical normalisation condition is 



' dz 



(34) 



4. Vertex residues 

In our subsequent analysis of deconfinement and chi- 
ral symmetry restoration we will display the response to 
changes in temperature of the residues connected with 
the TT- and cr-meson screening masses in, respectively, the 
pseudovector and pseudoscalar vertices, and the scalar 
vertex. 

In the isospin-symmetry limit, which we employ herein, 
there is no cr-meson pole in the vector vertex. On the 
other hand, the pseudovector vertex does exhibit a pion 
pole and the residue is the pion's leptonic decay constant. 
At T 7^ the expression for the relevant decay constant 
may be derived following Refs. psi. ISTj: 

^ = ^^[S^-2F^]^?G?£, (35) 

where "Kfe given in Eq. ()15c|) . This is actually . The 
expression for /J is different and of less interest herein. 

The pseudoscalar vertex also exhibits a pion pole. Its 
residue is [H, [H^j 



— [Et, IQe + Ft ^f] 5 



(36) 
The 



where the kernels are given in Eos. (|15ap . 
product r^/^ defines the in-pion condensate 

The scalar vertex exhibits a a-meson pole. Its residue 
may be derived following Refs. [stI. [soj: 



47r2 



cr I 1 



(37) 



where is given in Eq. (pS)) . This residue can be used 
to express the in-cr- meson condensate [59| . 
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C. Diquark correlations 

The relevance of rainbow-ladder truncation mesou- 
BSEs to baryon Faddeev equations is explained, e.g., 
in Sect. 2.1 of Ref. [2l|; namely, in this truncation one 
may obtain the mass and Bethe-Salpeter amplitude for a 
colour-antitriplet quark-quark correlation (diquark) with 
spin-parity from the equation for a J~^-mcson in 
which the only change is a halving of the interaction 
strength [l9| . The flipping of the sign in parity occurs 
because intrinsic parity is opposite for fermions and an- 
tifermions. N.B. Only scalar and axial-vector diquark 
correlations are needed for the ground-state nucleon and 
A because these correlations have the same parity as 
those baryons and masses which are lower. 

We remark that the rainbow-ladder truncation gener- 
ates asymptotic diquark states. Such states are not ob- 
served and their appearance is an artefact of the trun- 
cation. Higher-order terms in the quark-quark scat- 
tering kernel, whose analogue in the quark-antiquark 
channel do not materially affect the properties of vec- 
tor and flavour non-singlct pseudoscalar mesons, ensure 
that QCD's quark-quark scattering matrix does not ex- 
hibit singularities which correspond to asymptotic di- 
quark states [3^ [G^I . Studies with kernels that ex- 
clude diquark bound states nevertheless support a phys- 
ical interpretation of the masses, , obtained us- 
ing the rainbow-ladder truncation; viz., the quantity 
^{ggV^ '~ '^l'^iii)jp interpreted as a range over 
which the diquark correlation can propagate before losing 
its identity through fragmentation. 

Following these observations, it is straightforward to 
infer the BSEs for diquark correlations from the formu- 
lae in Sec. Ill Bl The Bethe-Salpeter amplitude for a spin- 
parity diquark is equivalent in form to that for a 
meson; e.g., the 0+ diquark is described by an amplitude 
ro+(Qo) and 



r[r+(Qo) 



ro+(Qo)ct 



(38a) 



^757-Qoi^o+(Qo), (38b) 

where C = 7274 is the charge-conjugation matrix [see 
Eqs.(|D10|, ([UTT]) ]. The amplitude V^+{Qo) and the 
mass of this correlation are obtained from an equation 
with the same form as that for the pion except for inclu- 
sion in the kernel of a multiplicative factor of 1/2. (See, 
e.g., Eq. (23) in Ref. [13.) 

As we are concerned with just the nucleon and A, only 
one more BSE is needed; namely, that for the axial- vector 
diquark. The Bethe-Salpeter amplitude for this correla- 
tion is constructed from 



rf+(go) 



74S^+(Qo) 

7_i£'i+ (Qo) 



(39) 



The mass of the longitudinal component of the correla- 
tion is the solution of 
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FIG. 1. Poincare covariant Faddeev equation, the basis for our 
computation of baryon screening masses. is the Faddeev 
amplitude for a baryon of total momentum P = Pq + Pd- The 
shaded region demarcates the kernel of the Faddeev equation, 
Sec.jB] in which: the single line denotes the dressed-quark 
propagator. Sec. Ill Al F is the diquark Bethe-Salpeter ampli- 
tude, Sec. lII Cl and the double line is the diquark propagator, 
Eqs. ([Bll), (fCej) . 



where is given in Eq. ((T9| ; and that of the transverse 
component from the same equation except that one omits 
in mapping — > !}(f . The amplitudes are canon- 
ically normalised as follows: 



1 -'.2d 



w-r-'^ii"^' (z) 



(41a) 
(41b) 



(40) 



D. Baryon Faddeev Equations 

We base our description of the dressed-quark-core of 
the nucleon and A resonance on the Faddeev equation 
introduced in Ref. d^], depicted in Fig.[TJ and since stud- 
ied extensively at T = (see, e.g., Refs. HH, [H [H IH- 
|68|). The phrase "dressed-quark-core" means that, con- 
sistent with the rainbow-ladder treatment of mesons and 
diquark correlations, we deliberately omit contributions 
to baryon masses that arise from resonant (meson cloud) 
contributions. The nature and implications of this omis- 
sion at T = are detailed, e.g., in Ref. [2l| (particularly. 
Sec. 4.2) and Ref. [13 (particularly. Sec. 3.3). 

With the loss of 0(4) invariance at nonzero temper- 
ature, a complete description of J = i, | baryons be- 
comes complicated. In the case of the nucleon, as with 
the dressed-quark in Eq. ^ , the complexity begins with 
a new structure in the propagator; and then one must 
account for the separation of the axial-vector diquark 
into two components. Notwithstanding this, a care- 
ful symmetry-preserving formulation using the contact- 
interaction could conceivably yield a tractable albeit clut- 
tered problem. For the A-resonance, on the other hand, 
the complexity becomes extreme [69|. 

In this first study, therefore, we choose to work with 
the zero temperature Faddeev equations modified in a 
simple manner. Namely, in deriving the Faddeev equa- 
tions we: ignore T-dcpendcnt modifications of the prop- 
agators and amplitudes; and then realise T-dependcnce 
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in the resulting equations by replacing the d/4 relative- 
momentum integral by a Matsubara sum, employing the 
T-dependent mass of the dressed-quark and the screen- 
ing masses of the diquarks (transverse mode, in the case 
of the axial- vector) , and the appropriately matched T- 
dependent diquark Bethe-Salpeter amplitudes. The pro- 
cedure is exemplified via the A in App.|B]and the nucleon 
Faddeev equation is described in App.lCl 



III. CONFINEMENT LENGTH AT NONZERO 
TEMPERATURE 

Recall now the rightmost expression in Eq. ([5]) : a finite 
value of Tir = l/Air « 0.8 fm implements confinement by 
ensuring the absence of quark production thresholds in 
all processes We expect that QCD exhibits dc- 

confinement at some T = > 0, whereat the production 
thresholds reappear, as illustrated in Ref. ^M,- Here we 
therefore introduce a dynamical mechanism that makes 
Til- temperature-dependent. 

Whilst more sophistication is required in general [70| . 
if one works in the rainbow-ladder truncation, then chiral 
symmetry restoration and deconfinement may be studied 
in the chiral limit by using the auxiliary-field effective ac- 
tion [zllzl, which we'U denote by J?. At T = 0, realistic 
models of QCD's gap equation support a DCSB (Nambu- 
mode) solution and a chirally-symmetric (Wigncr-modc) 
solution^ The difference 



'Bwn{T) := Jl[Wigner] - .^^[Nambu] 



(42) 



measures the relative stability of these different modes 
0: when 'Bw^^ is positive, the Nambu mode is dynam- 
ically favoured. The difference in Eq. (|42)) evolves with 
T; and the temperature at which it vanishes defines the 
critical value, T = Tc, for chiral symmetry restoration. 
Chiral symmetry restoration and deconfinement are si- 
multaneous in extant DSE studies; viz., Td = (see, 
e.g., Refs. [3,|4l|4l[7i-[i3)- Hence, following an idea in 
Ref. B, we define 



nr(T) 



(43) 



which ensures that the confinement length-scale diverges 
when chiral symmetry is restored in the chiral limit. 

Our precise implementation of Eq. (j43|) is described in 
connection with Eq. (|53|) . It uses the fact that, with a 
dressed-quark propagator of the type in Eq. ([2]), the ex- 
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FIG. 2. Temperature dependence of the dressed-quark mass 
in the chiral limit. Solid curve - the standard (positive) 
Nambu solution, Af+; long dashed curve - negative Nambu 
solution, which necessarily partners the positive solution; dot- 
ted curve - Wigner-mode solution; and short dashed curve - 
{-dM+/dT). 



plicit form of 'Bwn{T) is readily evaluated [li 



•BwHiT) = 2N,Nf / In 

JlAp I 

pf + mMN 



Pt 



■Pi 
Pf 



- mlVL 



w 



Pt 



Ml 



Pt 



(44) 



where pf = p'^ + ujf and Nf = 2, N^^ 3. Equation (gl]) 
possesses an ultraviolet divergence; and our regularisa- 
tion procedure is explained in connection with Eq. (jSOp . 



IV. RESULTS 

A. Phase transition 

In order to explore chiral symmetry restoration in our 
symmetry-preserving regularisation of the contact inter- 
action, we solved Eq. © in the chiral limit using the val- 
ues of aiR, Tir, Tuv specified in Table ID The solutions 
are depicted in Fig. [2] It is plain that chiral symmetry is 
restored via a second-order transition at 



T° = 0.215 GeV, 



(45) 



which is the temperature at which the chiral susceptibil- 
ity, X ■= —dM/dT, diverges in the Nambu phase. 

Denoting the positive and negative Nambu solutions 
by N±, respectively, then it is not surprising that in the 
chiral limit 



^ The full pattern of solutions to the gap equation is described in 
Refs. |73|477|1 . However, as explained in Ref. l77t . when discussing 
phase stability it is sufficient to consider only the simplest Nambu 
and Wigner solutions. 



VT : 'BwN^ (T) - 'BwN^ (T) = . (46) 

The value of T° in Eq. gS]) is between 20% and 40% 
too large when compared directly with that obtained 
in numerical simulations of lattice-regularised QCD [sTI - 
dll. It is notable, however, that r° = 0.234 to°, where 
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FIG. 3. Temperature dependence of the dressed-quark mass 
at m = 7MeV, the physical hght-quark current-mass in our 
model. Solid curve ~ the standard (positive) Nambu solution, 
M+ ; long dashed curve - negative Nambu solution; and dotted 
curve - Wigner-mode solution. 



TOp = 0.919 GcV is the model's zero-temperature chiral- 
limit value for the p-meson mass. Measured in these 
units, our value of T° is between 0% and 15% too large, a 
window that is typical of the rainbow-ladder truncation 
[ll,[2i[3l|. 

In Fig. [3] we depict the T-dependence of the simplest 
solutions of the gap equation at the physical value of the 
light-quark current-mass. By contrasting Figs. [5] and |31 
it becomes evident that the chiral symmetry restoring 
transition is replaced by a crossover at nonzero current- 
mass. In order to implement Eq. (|43|) . we must ask how 
then to assign a unique critical temperature for light- 
quarks with m 7^ 0? 

An answer is suggested by the behaviour of the N_ 
and W solutions in the Figure; and explained via a thor- 
ough consideration of the nature of the gap equation's 



solutions, as described in Ref. |77| . Expressed simply, in 
the neighbourhood of m = 0, DCSB is manifested in the 
simultaneous existence of three solutions to the gap equa- 
tion. When T reaches a value such that just one solution 
remains, explicit chiral symmetry breaking has come to 
dominate in the solution of the gap equation. For m > 0, 
therefore, is defined as the merging temperature of 
the and W solutions, which may readily be located. 
It is the common temperature at which the chiral sus- 
ceptibility diverges when evaluated with the W and 



phases or, equivalently, the solution set 

{rr} = {r>o|®wjv_(T) = o}. 



(47) 



In order to define this set, the integral must be regu- 
larised. We accomplish this by first noting that 



5m 



"SwN-iT) = {qq)w - {qq)N. , 



where 



{qq)p = -m,Nf [ 

Jlc 



Mp 



Idp Pi 



Ml 



(48) 



(49) 
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FIG. 4. Upper panel - Temperature dependence of the 
energy-density difference 'BwN^iT) at m = 7MeV. Lower 
panel - Fourth-root of that quantity: as noted elsewhere Q, 
"^vvN = 0) = 0.17 GeV is commensurate with the energy- 
difference assumed in bag-like models of baryons. 



an expression which can be regularised via Eq. (P. Now, 
at each temperature there is a current-quark mass, TOc, 
such that the W and 7V_ solutions merge and hence 
{'iq)\v ~ {qqYNl ■ Thus the regularised energy-density 
difference at current-mass m may be reconstructed as 



•BwNiT) 



dt 



{qq)w^' - {qq)N=' 



(50) 



At m = 7MeV, as illustrated in Fig.lH Eq. (j47l) yields 

(51) 



T^" 0.197 GeV = 0.212 < r° . 



The behaviour illustrated in the Figure does not de- 
pend sensitively on the interaction chosen: it is typical of 
a second-order symmetry-restoring transition (compare, 
e.g., the results in Refs. jsl. lil [sl - lsTt ) . 

We note that, evaluated with the 7V+ solution, the chi- 
ral susceptibility is maximal at 



rpn 



0.221 GeV 



0.238 = 1.12 



"X —^■--^^■■-p — ^.i^^c . (52) 

The preceding discussion clarifies Eq. and enables 
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FIG. 5. Calculated T-dependence of several quantities that 
are commonly used to illustrate the evolution of hadron 
properties through the chiral symmetry restoring transition: 
Solid curve - dressed-quark mass; dotted curve - pseudovec- 
tor component of the pion's Bethe-Salpeter amplitude, F-^, 
Eq. ((9|; and dashed curve - pion's leptonic decay constant 
fn, Eq.lSll). 



FIG. 6. T-dependence of the residues of the n- and a-mesons 
in, respectively, the pseudoscalar and scalar vertices: solid 
curve - r-ir in Eq. p6]): and dashed curve - ra in Eq. ((37]). 



haviour of in the neighbourhood of T"' is an artefact 
arising from interference between gso{T) and Tii.(r). 



a concrete implementation; viz., for light quarks: 



1/4 

"WN- 



(0) 



1/4 



(T) 



(53) 



All results described below are obtained by using this 
T-dcpendent infrared length-scale in the computation of 
the functions C\ Cf, . 

In Fig. [5] we depict the temperature dependence of 
quantities that may all be considered as equivalent chiral 
order parameters. The novelty, perhaps, is which is 
the pseudovcctor component of the pion's Bethe-Salpeter 
amplitude, Eq. (j9l). As remarked above: a pseudoscalar 
meson must possess components in its Bethe-Salpeter 
amplitude that may be described as pseudovector in char- 
acter ; and these pieces materially influence the T = 
physics of pseudoscalar mesons [2^, IH, H^, [3l|, . With 
increasing temperature, however, the strength of these 
components diminishes until, at in the chiral limit, 
they disappear [s^]- It is that outcome which forces /jr 
to vanish. This is abundantly clear when analysing the 
contact interaction: consider Egs. p^ . and recog- 
nise that the expression for , Eq. ([55)) , is proportional 
to the equation for F.^. That is equivalent to M as 
an order parameter is also plain: I^ef cx M [see 

Eq. (jl5bl) ] and hence the driving term for F^^ vanishes 
with M. 

In Fig. [6] we report the temperature dependence of the 
residues connected with the tt- and cr-meson poles in, 
respectively, the inhomogeneous pseudoscalar and scalar 
vertices. Qualitatively, the behaviour is similar to that 
depicted in Fig. 6 of Ref . jSTj : chiral symmetry is re- 
stored. Herein, however, the T = magnitude of the 
splitting Tc — is greater because we have modified the 
cr-meson BSE by the inclusion of gso in Eq. (^51) . The bc- 



B. Screening masses: mesons and diquarks 

We are now in a position to compute and report the 
T-dcpendcnt screening masses for the two-valence-body 
systems whose Bethe-Salpeter equations are detailed in 
Secs.HlBl InCl 



1. Screening in J = channels 

In Fig. [7] we plot the T-evolution of the screening 
masses associated with J ~ systems. At T = 0, our 
implementation of the contact interaction produces the 
following bound-state inertial masses j27| : 



71" cr [qq]„+ [qq]„- 
mass (GeV) 0.14 1.29 0.78 1.37 



(54) 



Plainly, DCSB is expressed strongly in this part of the 
spectrum through a large splitting between parity part- 
ners. The splitting persists in the screening masses. 
They are approximately insensitive to temperature until 
T = r™ = Td, beyond which value the screening masses 
of the parity-partner correlations evolve rapidly toward 
near degeneracy: apart from current-mass effects, they 
may be called equal for T > 1.3 T™. 

Since our confinement mass-scale vanishes for T > 
T™ Td; i.e., Ai,.(T > T™) = 0, it is noteworthy that the 
correlations persist on this domain. Thus, deconfinement 
is not expressed in the absence of strong correlations in 
these channels. This is because fermions at nonzero tem- 
perature are characterised by an additional mass-scale: 
ujQ = ttT; and so long as 2[ujq+M{T)] exceeds the correla- 
tion's screening mass, the correlation will persist. These 
features were exposed in Ref. [STj . 
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Considering the construction described in connection 
with Eq. ^ , one should expect a dramatic expression 
in the real-time propagator of the discontinuous deriva- 
tive exhibited at by the lightest chiral-limit screening 
masses in the tt and a channels (see Fig. 3 in Ref. [s^] ) , 
and its smoother remnant at m > 0, which is evident in 

Fig.m 



2. Bound states? 

Given these observations, we judge that the behaviour 
of the screening masses does not necessarily mean that 
bound-states continue to exist in any given channel at 
T>T,. 

There are illustrative examples; e.g., Ref. [ssf . which 
considers the scalar and pseudoscalar channels. In ideal 
rainbow-ladder truncation (i.e., with 530 = 1): the iner- 
tial mass of the scalar-meson is — 2M in the chi- 
ral limit, and the pseudoscalar mass is zero. As the 
T ~ interaction strength is reduced to some criti- 
cal value (e.g., aiR is reduced), the dressed-quark mass 
M — >• and hence 0, thus becoming degenerate 

with TTiTT ==0. As the interaction strength is reduced 
still further, the bound-states disappear. It is possi- 
ble, therefore, that when some mechanism suppresses the 
interaction-strength to a sufficiently large extent, dynam- 
ical mass generation is impossible and no true bound- 
states are supported. 

With these observations in mind, we conjecture that, 
when they may reasonably be defined, the inertial masses 
of all hadron bound-states are proportional to M{T) for 
T < Td] and that no bound-states persist above Td- 
This suggestion is less surprising if one states it thus: 
as the critical temperature characterising a second-order 
phase boundary is approached from below, all correla- 
tion lengths diverge and real-time correlation functions 



acquire power-law behaviour. In our case, the primary 
correlation length is ^(T) ~ ru.(T) ~ l/M{T), the di- 
vergence of which forces all related correlation-lengths - 
the inverse of hadron inertial-masses - to diverge as well. 
Notably, for states with a significant hadronic width, a 
vanishing mass might be very difficult to distinguish em- 
pirically from marked spectral broadening: M — > rii — ioj, 
with oj/m ^ 1. We note in connection with the lat- 
ter that since rainbow-ladder truncation omits resonant 
contributions to bound-state kernels, it is likely to be- 
come a quantitatively inaccurate approximation if used 
to compute the dynamical evolution of bound-states in a 
medium with copious numbers of thermal pions. 

The possibility that hadron inertial-masses decrease in- 
medium has long excited interest 18^ ; and many analyses 
have explored this possibility [Tsl-IOOl l9l| . At present the 
notion is empirically neither confirmed nor invalidated, 
although an explanation of data does not require this ef- 
fect in its simpler forms [92[ . A reduction in hadron iner- 
tial masses in the neighbourhood of Tc is consistent with 
results from some applications of sum rules but spectral 
broadening via hadron-hadron interactions in-medium is 
also very important Much remains to be learnt 

in this area, at least in connection with the role and man- 
ifestation of gluon-quark dynamics. 



3. Screening in J = 1 channels 

We plot the T-evolution of screening masses associated 
with J ~ 1 systems in Fig. [51 At T = 0, our implemen- 
tation of the contact interaction produces [27| : 



P ai {qq}i+ {qq}i- 
mass (GeV) 0.93 1.38 1.06 1.45 



(55) 



For T > 0.3 T™ a separation is apparent between 
the screening masses of the transverse and longitudi- 
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FIG. 9. Screening masses of the nucleon and A dressed-quark 
cores: solid curve - nucleon; and dashed curve - A. 



FIG. 10. Difference between screening masses: soUd curve - 
A and nucleon dressed-quark cores; and dashed curve - axial- 
vector and vector diquark correlations. 



nal modes. The subsequent behaviour of the screening 
masses for the transverse modes follows the pattern set 
by J = systems. They are weakly sensitive to tem- 
perature within the confinement domain, a result found 
previously in the algebraic model of Ref. [Oal (see Eq. (21) 
therein) ; and chiral symmetry restoration is again evident 
for T > 1.3 T™, apart from current-mass effects. 

The longitudinal modes behave differently, however. 
Their screening masses all increase markedly with tem- 
perature, a result readily understandable from Eq.(22) 
in Ref. [96[, such that, with the exception of rjipW , they 
are greater than 2uja at the dcconfincmcnt temperature, 
Td = in Eq. (fSTj) . In these three channels, owing to 
the additional repulsion produced by the ^'"(T) term in 
the Bethe-Salpeter kernels [see, e.g., Eqs. (|32|) . (|33| ]. this 
is sufficient to dissolve the correlations. 



C. Screening masses: nucleon and A 

Having determined the behaviour of the screening 
masses and Bethe-Salpeter amplitudes of the diquark 
correlations, we are now positioned to compute and re- 
port the temperature dependence of the dressed-quark 
cores of the nucleon and A as they are described by the 
Faddeev equations detailed in Sec. Ill Dl The T = in- 
ertial masses are listed and discussed, respectively, in 
connection with Eqs. (|C19|) . (|B2ip . Their evolution into 
temperature-dependent screening masses is depicted in 
Fig.m correlations persist in both channels for all val- 
ues of T; and for T > T™, the screening masses increase 
linearly with temperature but always lie below SttT. 

It is interesting to analyse the splitting between the 
A and nucleon screening masses. In the T = study 
of Ref. [Ill it was shown that, as a function of current- 
mass, mA — m^r is linearly proportional to the split- 
ting — '^[gg]o+ • Figure [TU] demonstrates that 
a similar correspondence holds at fixed current-mass 
with increasing temperature: given the T-dependence of 



™{'3'3}i+ ""^[9910+ ' ^^^^ that of m./y^—mN is approximately 
the same. 

The existence of diquark correlations within baryons is 
a dynamical outcome of the strong interaction between 
quarks. Whether one exploits this feature in order to 
develop an approximation to the quark-quark scattering 
matrix, as we do herein, or chooses instead to eschew the 
simplification it offers, the outcome is the same It 
follows that axial-vector diquark correlations are domi- 
nant within the A. The nucleon, on the other hand, pos- 
sesses both scalar and axial- vector diquarks; and the nu- 
cleon's Faddeev amplitude expresses the relative strength 
of these different correlations within the nucleon. The 
T = result is presented in Eq. (|C20|) : the scalar diquark 
is found with 72% probability. (The significance of this 
result for the hadron spectrum is described in Ref. (27[.) 

It is natural to consider how the ratio evolves with tem- 
perature in the correlations that persist above Td = T™. 
This is depicted in Fig.[TT] whilst there is a perceptible 
evolution in the apportionment of strength between the 
two axial-vector diquark components, with the relative 
probability switching at T = Td, the probability of find- 
ing a scalar diquark correlation is almost temperature 
independent. This is a curious outcome, which owes to a 
relative similarity in the evolution of both Eq+ , Ei+ and 
TOo+ J 1 ^-s illustrated by Fig.[TUl In contrast. Fig. [5] 
shows that the structures defining the pscudoscalar corre- 
lation, and hence the scalar diquark, both change signif- 
icantly above Td- It will be interesting to learn whether 
the behaviour in Fig. [11] survives a more sophisticated 
treatment of the T Faddeev equation. 

One might also ask about parity partners in the baryon 
sector. Based upon the expressions in Ref. [l^l and their 
relationship to those derived herein, we anticipate that 
chiral symmetry restoration will lead to degeneracy be- 
tween = ^^,5 screening masses and, separately, 
J"^ = 1^, I masses. 
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FIG. 11. Evolution of the nucleoli's Faddeev amplitude with 
temperature: solid curve - scalar diquark component; and 
dashed, dotted curves - the two distinct axial-vector diquark 
structures. [See Eqs. ((Cl]), ((C2|l .] 



FIG. 12. Temperature evolution of the dressed-quark axial- 
charges in Eqs. (|64|l . evaluated with m = 7MeV: solid curve 
- g\^; and dashed curve - (j^^. 



V. GOLDBERGER-TREIMAN RELATIONS 

In- vacuum, the axial- vector vertex, r5^(fc_|_,fc), 
k + Q, is the solution of 



167raiR f dH 



3ml J (27r)' 



7a5(t + 0)r5p(Q)5(t)7a. (56) 



It satisfies the axial-vector Ward-Green- Takahashi iden- 
tity, which reads, in the chiral limit: 

P^.r5^.ik+, k) = 5-^(fc+)z75 + ilzS-\k) . (57) 

Translating the general relations in Ref. [i^ to our 
symmetry-preserving formulation of the contact interac- 
tion, Eq. ((57|) entails the following chiral-limit quark-level 
Goldberger-Treiman relations: 



7^0 _ 71 rO _ 77O 



(58) 
(59) 



where, as above, the superscript "0" denotes a quantity 
computed in the chiral limit and we have used the fact 
that, in this limit and in the neighbourhood of = 0, 
the axial-vector vertex has the general form: 



r«^(fc+, fc) = 757p^^,(Q) + ^V'.niQ) . (60) 



Here, r"(Q) is the canonically normalised pion Bethe- 
Salpeter amplitude and we have defined a dressed-quark 
axial-charge HH 



9% = F'AiQ - 0) . 



(61) 



With the loss of 0(4) invariance at nonzero tempera- 
ture, the contact-interaction axial- vector vertex takes the 



general form (Qo = {Q, 0}) 



r5/.((3o) = 75 



27||F/(go) + 2z^i?^(go), 
(62) 

where 7 =: 7_l + 7|| with Q • 7|| = Q • 7- In this case, 
the axial-vector Ward-Green- Takahashi identity does not 
place a tight constraint on F^^{Qo). However, the exis- 
tence at T 7^ of a pscudoscalar correlation with zero 
screening mass entails: 



E 

M°F, 



-°(o) 
-°(o) 



M°F^"(0) + 2/>0(0) 



/Oi?0(0) = M", (63a) 
f^F^iO), (63b) 
A/° = /,X(0); (63c) 



and vicc-vcrsa. 

As indicated in Eq. (pT|) , the regular parts of the axial- 
vector vertex may be identified as axial-charges of a 
dressed-quark. The T-dependence of these charges is de- 
scribed by the following formulae: 



giJT) = F\0) 



1 



l + ^^-(0)' 
1 



(64a) 
(64b) 



with the kernels given in Eqs. ([50]) . I\'3'3\i : and displayed 
in Fig.[T21 The behaviour is consistent with the T ^ 
Goldberger-Treiman relations, discussed in association 
with Eqs. (|62p . ((63)) . The charges are identical and less- 
than one for wT < M(0); and also essentially indepen- 
dent of T on this domain (consistent with Ref. (9^, al- 
though the context therein is different). However, 

and g\ become distinct on ttT > A/(0). The transverse 
charge remains below one until T « T™ whereafter it 
rapidly approaches unity, as dictated by Eq. (|63cp and 
the restoration of chiral symmetry. Like the screening 
masses, however, the longitudinal axial-charge is far more 
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sensitive to temperature than g^^. At T = T™, g'j^^ has 
dropped 20% from its T = value. With increasing tem- 
perature thereafter, it approaches unity from below; but 
only slowly: e.g., even at T = 2T™ it has not returned 
to its T = value. 

The connection in-vacuum between the drcsscd- 
quark's axial-charge and that of the nuclcon is discussed 
in Ref. [3l| . We revisit that here in the context of nonzero 
temperature. In this case the free-field nucleon spinor is 
defined via 

= [i^.p + i 74w„ +Mn]u„{P), (65a) 

= S„ {P)[lJ-P + i-fiiOn + Mn] , (65b) 

where Mjv is the fermion's inertial mass. It follows that 

Un{P')-/5^-{P'-P) UmiP) 

= u„(f'')75[-74('^« - t^m) - 2iMN]Um{P) ■ (66) 

Now consider the T 7^ extension of the nucleon's 
axial- vector current with momentum Qq entering the ver- 
tex: 

If one specialises to the chiral limit, then, analogous to 
the dressed-quark case, the existence of pseudoscalar cor- 
relation with zero screening mass entails 

UniP')l59piQ^)UniP) 

'^=' 2^g°^^(0) S„(P') ^75U„(P) , (68) 

where g^j^j^{0) is a normalisation factor that ex- 
presses the = value of the seven-point function 
Un{P') ij5Un{P). As a consequence of the nuclcon-lcvel 
axial- vector Ward-Green- Takahashi identity, one has 



= iQ-fiiQ) (69) 
= u^{P')ii^[i-Qgi\Q^) + Q^gUQ^)]u.n{P)m 



<3 ~o , 



2un{P') 75[Mjv.9i°(0) - flg%^]un{P) ; (71) 

and hence, in the chiral limit, at all values of tempera- 
ture: 



M]^37(0)-/X^A.(0). 



(72) 



In the chiral limit, = at T°. As with r^(Q), 
VT the normalised value of ?2„(P') i^5U„{P) is finite at 
= 0; i.e., < ,9°^jy(0) < 00. Consequently, 



lim M^g;i°(0)=0. 



(73) 



Following the reasoning in Ref. [31| , since g\^ is always 
nonzero and, indeed, si°(T = T°) « 1, Eq. ([731) is not 
achieved by changes at the level of the dressed-quark- 
axial- vector vertex. Hence a vanishing of g^{T^) would 
require extraordinary and precise cancellations amongst 



the terms that constitute the axial-charge matrix element 
associated with the correlation in the nucleon channel; 
namely, between the various contributions arising from 
the angular momentum correlations within the Faddeev 
amplitude. Owing to the power of symmetries in quan- 
tum field theory, this is conceivable. However, it is un- 
likely given the behaviour in Fig. 1111 i.e., that the am- 
plitude describing the T 7^ correlation in the nucleon 
channel is only weakly sensitive to T . It is more prob- 
able, therefore, as argued in Ref. [U , that Eq. ([75)) is 
achieved via dissolution of the nucleon bound-state at a 
point of coincident chiral symmetry restoration and de- 
confinement, with 



lim M?, 

T^(T°)' 



0, 



(74) 



so that, beyond Tc, g^^ represents just the normalisa- 
tion of a seven-point function that is associated with a 
strong screening-correlation but not a bound-state. The 
discussion in Sec. lIV B"2] anticipates Eq. ([71)) . 

The above discussion is dubious in the neighbour- 
hood of if the inertial mass of the nucleon acquires 
a large imaginary part within this domain; i.e., if M]^ — 
mPjq — iw^. with Lo'j^/rnP]^ > 1. In this case, Eqs. ([55)1 
become poor approximations. In fact, there is no sense 
in which one may employ notions of an asymptotic nu- 
cleon state and the definition of each of the vertices em- 
ployed above must be revised significantly. But this is 
also the content of Ref. [3l|: under these conditions, the 
Goldbcrgcr-Treiman relation is made moot by bound- 
state dissolution. 



VI. SUMMARY AND PERSPECTIVE 

Working at leading-order in a symmetry-preserving 
truncation scheme for the Dyson-Schwinger equations, 
we extended a confining formulation of a vector x vector 
contact-interaction to nonzero temperature. This frame- 
work proved useful in the study of a wide range of phe- 
nomena at T 0; including hadron masses and form fac- 
tors. We therefore expect that, interpreted judiciously, 
results obtained at T > should represent a fair guide 
to related hadron properties on this new domain. 

In formulating the interaction at nonzero temperature, 
our treatment of the gap and Bethe-Salpeter equations is 
standard. However, our formulation of the baryon Fad- 
deev equations is novel. Although, in common with pre- 
vious continuum studies, our treatment fails to express 
the full complexity of J = i, | states at T > 0, it does 
improve upon preceding analyses in a number of ways. 
For example, by: including axial- vector diquarks in addi- 
tion to the scalar correlations; improving materially upon 
the implementation of the widely-used static approxima- 
tion for the Faddeev kernels; and allowing dynamically 
for the expression of chiral symmetry restoration and, 
importantly, deconfinement. 
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Turning to the results, via the gap equation wc found 
a second-order chiral-symmetry restoring transition at 
T ~ — 0.23 mp in the chiral hmit, which becomes 
a cross-over at nonzero current-mass. Notwithstanding 
this change at m ^ 0, by capitahsing on the modern 
understanding of the gap equation and the nature of its 
solutions, we were still able to define a single tempera- 
ture, T™ = 0.21 TTip, whereat dynamical chiral symmetry 
breaking (DCSB) is no longer effective. 

We thereafter explained and implemented a dynamical 
mechanism that ensured deconfinement at T4 = T™. It 
is distinguished by the feature that, whilst no kernel or 
process exhibits production thresholds for coloured states 
when T < T^, all do for larger temperatures. Despite 
this, in all T = bound-state channels, strong corre- 
lations persist for T > Td- Furthermore, in the mass 
spectrum defined by the associated screening masses, de- 
generacy between parity partners is apparent for T > 

Since we retained axial-vector diquark correlations, we 
were able to simultaneously study nucleon and A-baryon 
properties. The splitting between screening masses in 
these channels evolves with temperature in a manner that 
is approximately proportional to the splitting between 
the axial-vector and scalar diquark masses. Curiously, 
we found that the scalar-diquark content of the corre- 
lation in the nucleon channel is almost independent of 
temperature: it is 72% on T < 2Tc, which is the highest 
temperature considered herein. 

Notwithstanding our results for screening masses 
and the associated correlation amplitudes, we argued 
that there are reasons (amongst them, the nucleon's 
Goldberger-Treiman relation) to suspect that, when they 
may reasonably be defined, bound-state inertial masses 
vanish as T — >■ T^; and that this is a signal of bound- 
state dissolution at the deconfinement temperature. 

In closing we reflect briefly on the question of the 
veracity of our results. Our formulation is internally 
consistent and more systematic than preceding contin- 
uum studies of the properties addressed herein. We have 
furthermore eliminated two parameters used in earlier 
formulations of the contact-interaction baryon Faddeev 
equations at T = and left the other three untouched. 
The material weakness is that the rainbow-ladder trun- 
cation omits contributions to the gap and bound-state 
kernels which might fairly be described as meson-loop 
corrections. However, whilst such corrections are nec- 
essary in order to reliably determine critical exponents 
associate d wi th the transitions (see Ref. [3], p. S49, and 
Refs. [99l . ll00l |). we do not expect them to have a material 
effect on the quantities we've focused upon herein. This 
assumption can be checked. For the present, however, we 
note that our predictions, e.g., regarding the persistence 
of correlations above Tc and the nature of chiral sym- 
metry restoration and deconfinement, are broadly con- 
sistent with analyses of contemporary lattice simulations 
in those instances where a sensible comparison is possi- 
ble. It is interesting now to turn this analysis toward the 



questions we've raised herein; e.g., developing a contin- 
uum connection between the screening masses and the 
real-time structure of spectral functions. 
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Appendix A: Contact interaction 

The key elements in our analysis are the dressed-quark 
propagator, and the meson and diquark Bethe-Salpeter 
amplitudes. All are completely determined once the 
quark-quark interaction kernel is specified. We use 



(Al) 



where mc = 0.8 GeV is a gluon mass-scale typical of the 
one-loop renormalisation-group- improved interaction de- 
tailed in Ref. |lOl| . and the fitted parameter aiR = 0.937r 
is commensurate with contemporary estimates of the 
zero- mom entum value of a running- coupling in QCD 
10M07| . We embed Eq. ^ m a rainbow-ladder trun- 
cation of the DSEs. This means 



(A2) 



in the gap equation and in the subsequent construction 
of the Bethe-Salpeter kernels. 

Whilst the interaction in Eq. ()Aip may be viewed as be- 
ing inspired by models of the Nambu-Jona-Lasinio type 
0, our treatment is atypical. Used to build a rainbow- 
ladder truncation of the DSEs, Eqs. (|A1|) . (|A2p produce 
results for low-momentum-transfcr observables that are 
directly comparable with those produced by more sophis- 
ticated interactions, as illustrated in Refs. [2TI. [23l - [28j . 

In Table HI for reference, we report T = values of 
u- and s-quark properties, computed from Eq. ([6]). The 
ratio nis/ni, where to = ( niu + md)l'2-, is consistent with 
contemporary estimates [l08l | . The result Ms — TOs w 
Mq is typical jl09l Ill0l | and indicates that the additive 
impact of DCSB is nearly as great for the s-quark as 
it is for M, d- quarks. In general, however, Mf — to/ is a 
monotonically decreasing func tion of to/ , bounded below 
by zero as to/ — >■ 00 jllOl llll| . 

Wc have simplified the form of Eq. © by introducing 
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TABLE I. Dressed-quark properties, computed from the gap 
equation and required as input for the Bethe-Salpeter and 
Faddeev equations, and computed values for in-hadron con- 
densates [601 ~ at T = 0. The results are obtained with 
am = 0.937r and (in GeV) Ai^ = 0.24 , Auv = 0.905. (These 
parameters take the values determined in the spectrum cal- 
culation of Ref. [2ll|. which produces rUp = 0.928 GeV; we 
assume isospin symmetry throughout; and all dimensioned 
quantities are listed in GeV.) 

iriu rris nis/mu Mp Mu Ms Ms/Mu i^l^^ nl^^ 
0.007 0.17 24.3 0.36 0.37 0.53 1.43 0.241 0.243 0.246 



elucidated below.) The amplitude is obtained from the 
following Faddeev equation: 

/'^(4;g)<(g) 



iHf.(^o,^;Q)r(^;QX(Q), (B2) 



with A'o = -if) + Q, K = -£ + Q, = -m|^ and 



where 



the function 

T) = ST ^ / dr e-^^'^+'^''> 



1 — — 00 



dq e 



l. — -~-oo 



(A3) 



'dTe-^'2TMe-^^^^^')^, (A4) 



where 1^2 (a;) is a Jacobi theta- function jll2j |. It is 
straightforward to make the connection with zero tem- 
perature results once one appreciates that 



(A5) 



In deriving Bethe-Salpeter equations for vertices 
and bound-state amplitudes, the Ward-Green- Takahashi 
identities are crucial. At T 7^ 0, the identity in Eq. 
is necessary and sufficient to ensure they are satisfied, 
where 



3^-(^;T) 



d 1 1 
dr 2t 



2^2(0-"^"'^'). (A6) 



Using Eq. (|A5|) . it is straightforward to show that 
!^'"(^;T^0) =0. 



Appendix B: Faddeev equations for A-baryon 

Here we explain the origin of our simple Faddeev equa- 
tions for the baryons' dressed-quark-cores. 

Using a symmetry-preserving treatment of the contact 
interaction at zero temperature, the A-baryon Faddeev 
amplitude can be written 



(Bl) 



whe re u^ {Q) is a Rarita-Schwinger spinor, defined in 
Eq. (jD12[) of App.IdI (The reason for this simplicity is 



KtiK) = Zr^'ll , T,M<) = ^P. + ^ , (B4) 

"'1+ 



is the axial- vector diquark's propagator and Tj^{K) is 
its Bethe-Salpeter amplitude: 



(B5) 



At this point, one post-multiplies by up{Q] r) and sums 
over the polarisation index to obtain [Eq. (|D14p ]. 

k+{Q)R^p{Q) = A I -0^ ^U^^(eo,e;Q) A+{Q)R,f,{Q) , 

(B6) 

which, after contracting with (5^^, yields 



1 = 2tr 



d'^i 



OF^ f f 7pH7-(4-A')+M]7^ 
[-ij -i + M] [-17 -Q + ttia] 



£2 



2mA 



'-Tp.iK)R,^{Q) .(B7) 



(N.B. Here and below we suppress f{£o,Q), f{£,Q) on 
the left- and right-hand sides, respectively, because, sub- 
ject to our approximations, they will finally cancel.) 

Previous T — studies of Eq. (|B7P have used vari- 
ants of the so-called "static approximation," in which the 
quark exchanged in Fig.[Tl described by S''^(— £0 + K) in 
Eq. (jB3|) . is replaced by g\/M. This expedient is dis- 
cussed extensively in Sec. 4 of Ref. [2]| and Sec. 3.1 of 
Ref. [l^l . In starting with Eq. (|Bip , we assumed im- 
plicitly that a truncation of this sort would eventually 
be made, for only then does ijj^{l\Q) take such a sim- 
ple form. Indeed, in combination with diquark correla- 
tions generated by Eq. (jAl[) . whose Bethe-Salpeter am- 
plitudes are momentum-independent, the static approx- 
imation generates Faddeev equation kernels which are 
themselves momentum-independent and hence so are the 
Faddeev amplitudes. The consequent simplifications are 
the merit of the truncation. 

Unfortunately, the static approximation is inadequate 
at nonzero temperature because, in the chiral limit, the 
dressed-quark mass is expected to vanish for T > Tc- 
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In order to maintain context with the large body of work 
that has used the contact interaction pll. l23l - [28| . we must 
provide a reasonable alternative. 

To this end, consider that the right-hand-side of 
Eq. (HZl has the form 

dH N{£,£^ Kq,K,Q) 



(27r)4 [£2 + M2] [{i-Ko)^ + M2] [{-£ + Q)^ + m 



i+J 



(B8) 

where N{£, £~Kq, K, Q) is a numerator and we have used 
the relation (— £o + K) = (— £ + Kq). In proceeding from 
this point, one employs a two-variable (a, /?) Feynman 
paramctrisation to convert Eq. (|B8p into 



(B9) 



2 1 da dp a 
'o 

dH N{1 + TjpQ, l-{p- rip)Q, K, Q) 
(27r)4 [Z2 + ,s(A./2,m2^,a,/3,p,g2)]3 



where we have written Kq = pQ; i.e., explicated that 
the external diquark carries a fraction p of the baryon's 
momentum, defined i^p ~ a + aPp and also 

^B(M^m2+,a,/3,p,g2) = aM^ 

+ aml+ + [{a + a^p^) -{a + aPp)^]Q^. (BIO) 

Every term in the denominator is real, and hence it 
acts as a weight function whose maximum occurs when 
P = and 



iPp : 



'1 + 



2g2 



(Bll) 



We therefore define the integrand in the Faddeev equa- 
tion as the quantity obtained from Eq. (|B9I) with 



1 



- mj^ 



2Q2 



(B12) 



In so doing, we solve for the baryon Faddeev amplitude 
evaluated at a single value of the quark and diquark mo- 
menta; viz., respectively, = (1 ^ Pi+)Q^ Kq = pi+Q. 
(N.B. In a weak binding approximation, = 4M^, 

g2 = -m| = and hence pi+ = 2/3.) 

To be explicit, applying these rules in the present case 
maps Eq. (|B7p into the following Faddeev equation for 
the A-baryon: 

1 = / da dp a 
"lA 7o Jo 

dH 1 

2^,13 



xtro 



-M][-ij ■ Q + mA\Tpu{K)R,^{Q) 



(B13) 



The numerator in Eq. (|B13[) . and analogous numerators 
in the nuclcon Faddeev equation, produce many inner 



products. We treat them as at T = 0; viz., using the 
rules of Ref. [lH: 

{X-l){Y -l) ^PX-Y, X,Y = Q,K,Ko, (B14a) 

l-Q^ VQ\ l-K^ TjpQ\ I ■ Ko ^ ^?p/g^(B14b) 
K-P^pQ\ Ko-P^ PfQ^, (B14c) 



^0+ 



^0+ ' ^1+ 



(B14d) 



with p = 1 — r], r] — d + appf ; and pf = pi+ in Eq. (jB12p , 
if the external diquark is an axial- vector, or 



Pf 



•'0+ 



2g2 



=: Po+ 



(B15) 



if the external diquark is a scalar correlation. 

We now define the nonzero temperature Faddeev equa- 
tions by replacing the fi(74-integration in the expressions 
that arise by a Matsubara sum, so that 



/ 
/ 



dH 



dH P 



Cl{^;T), (B16) 
Pi'(^;T)-<,C2(^;T)], (B17) 



where 
tl{,;T) 



1 d_2 

2dc 



(B18) 

with d\^; T) defined in Eq. (gH). 

The procedure described above, applied to Eq. (jBlSp , 
yields the following Faddeev equation for the A: 

l = ^^(g' = -mi), (B19) 

with, using = ^AiM^ ,ml+ , a, P, pi+ , ~m\), 



^'^{Q' = -mi) 



= ^ j^dadPa |cr(qT) 
M + mA?7 



1+ 



6M + mA{2p + 8pi+ -677) 



mipl 



ruAi^P Pi+rj " ipl+i] 



"1+ 



2p,^p' ~3r^p')] 



2 m 4:p^pl+rim% 



'1+ 



(B20) 



Naturally, Ei+ = Ei+{T); i.e., the axial-vector diquark's 
Bethe-Salpeter amplitude is T-depcndent, as are M, 
mi+. 

At T = 0; using the definitions in Table HI this Faddeev 
equation produces 



mA(T = 0) = 1.42 GcV. 



(B21) 



17 



Equation (jB20|) represents an improvement over previous 
definitions of the static approximation because it elimi- 
nates one parameter yet produces the same mass (within 
2%) as that fitted using the extra parameter. 

We reiterate here that this equation describes the A- 
baryon's "dressed-quark-core" . It should and does, there- 
fore, produce a mass that lies above that quoted em- 
pirically for the A. Remarkably, indeed, the value pro- 
duced is almost identical to that inferred for the dressed- 
quark core via the dynamical coupled-channels analysis 
of Ref . |ll3j . Any theoretical framework that produces a 
stable (width zero) A-resonance with a mass near that 
quoted empirically is dubious. In this connection it is also 
noteworthy that m/^/mp ~ 1.53, a result which compares 
favourably with the experimental value: 1.59. 



Appendix C: Nucleon Faddeev equation 

We capitalise on Eq. (C.47) in Ref. [2l| to write the 
nucleon Faddeev equation as a 3 x 3 matrix eigenvalue 
problem for the scalar [ud] and axial- vector {ud} diquark 
correlations. To begin, the nucleon is described by a 
Faddeev amplitude: 



s{Q) 



uiQ), 



(CI) 



where 



siQ) = s(Q)Id , (C2a) 
flp(Q) = ai(Q)757A. + a2{Qh5Q,. , (C2b) 

with = —1 and u{Q) being the spinor introduced in 
Eq.jDl. 

The amplitude satisfies 

I.N 



d^i 



(27r) 



5W{^.}(4,^;g)V'M(Q)^(Q), (C3) 



with 



(C4) 



iWoo 3(5Woi), 
(5Wio)m -(^Wii)p^j 

where {l^ ^ Iq - K) 

= r"+(A-) S^{-i^) f "^(-/^o) S{i) A"^ {K) , (C5a) 
(%i). 

= v^;{K)s'^{-^^) f^\~K,) s{e) aIUk), (C5b) 

(5Wio),. 

= T"\k) S^{~£^) fl\-Ko) S{e) A°^ (K), (C5c) 
(^Vfii)^!^ 

= rl\K) S^{-1^) f i^(-Ko) S{1) AltiK) . (C5d) 



In Eas. (fC5a)) - (IC5c)) . 

a"+(a:) 



(C6) 



is the scalar diquark's propagator and (A') its Bethe- 
Salpeter amplitude: 

ro+ {K)C^ = z75i?o+ (K) + ^757 ' KFo+ (K) . (C7) 

At this point one can follow the procedure in Apd.IB] 
and thereby arrive at the following eigenvalue problem 
for Q2 = 



I'M- 



ai(Q) 
a2{Q) 



(C8) 



where the kernel is a 3 x 3 matrix: 



^(Q) 



Ass '^Asai 'JAsa2 

Vis Viai Aiia2 

_frll —"it'll 

-'\l2S -'\12a1 -'\L2a2 



(C9) 



The entries in Eq. (jC9p are described below. N.B. The 
Bethe-Salpcter amplitudes, dressed-quark mass and di- 
quark masses are T-dependent. 

AT™ = [eI+ %ee + -^0+ ^0+ 'Kef 

+fI^%ff\. (ClOa) 

Xee^ I dadfia{c"^(s-T)+t2(s\T)\{M 
Jo 

+mN [po+ - v] ) {M + mNr])-^]}, (ClOb) 
Kef = J^J/^d/3a {r,^{,;T){po^ - 2p) 

-r"(q T)[2(po+ + p){M + m.NvKM 
+mNipo+ - r])) + i{po+ - 2p)] } , (ClOc) 



Kff 



dadpap{ct{,;T)-C2{<^;m 
+ 2(M + mjv(po+ ' V)){M + m^yy)]} , (ClOd) 



where = ,ml+, a, P, p„+ , -m%), p = I - -q, -q = 

a + aPpQ+ . 



Kfa, = ^[E^^KEa.+F^+^a,], 
1 



= / dadpa 1 3^1 {^\T) + C2 (t;^) 

-I- Sm%ri{po+ - 77) + Pq+Mttin 
2p'^po+Mm% ' 



(Clla) 



sivr ~ 3<r 



(CUb) 
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'Kpa 



■mNPo+ 
2M 



y dadp a 1 Ci T) 


J- 1 9 







2A/2 + 6A/mArpo+ +^ 



TO 



1 + 



and 



01 



-^-j [i?o+ + Fa+ 'Kpai] 



(Cllc) 



(C12a) 

+ (Af - TOArr/)(TOAr(po+ - ??) - M)] } , (C12b) 



2M 



+2{M - mNv)iM + fUNiv - po+])] } , (C12c) 



where = ^B(A/2,TOj+,a,/3,po+, ''1+ 
m%p'^/ml+ , p=l-r],T] = a + a/3po+ ■ 



/•I 



(C13a) 



6to^+ 



+ C"(qr)(2A/2_2TO2,r^2_^)] 

+'^2"('?;r)pi + TOAr(A/ + 77mAr)| , (C13b) 
2to^+ + TO^p^+ _iu 



6to^^_ 



+4"(^;r)(2A//2-2m^r?2-,)] 
+ ^^2"('?; T)pi+mN{M + rjmN) 



(C13c) 
(C13d) 



where c = (;B{M^,ml+,a,P,pi+,-m%), p = 1 - 77, 77 
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-E 



[Eo+ 'Ka2E + -Fo+ ^02 -F 
1 r^2 /|^2»,2 

dad/3 a 



(C14a) 



,,-4to,,,,_.^^^^^ 



6to^+ 



1- 



-^G (^;r)] + -C2 (qT)(Af + 77TO,v) 

X [A/ + TOAr(577 - 3pi+) 



^02 F 



1+ 

^ f'dadPa^'-^^^^^^^^^^r^ 







(2A/ + tjitin)] } 

{ 



(C14b) 



+ 4"(^;T)(2Af2 _ ,)] + c;"(qT)[277TO^(3pi+ 
+r]) + 6MmNipi+ + 2r/) - '^'^^^1+ 



-1+ 



X (4772m^ + 677A/77ijv] | , (C14c) 

where <j = '^b{M'^ ,ml+, a, P, pi+ , -m%), p = 1 - 1], i] = 
a + aPpi+ . 



2ml+ 



'I. 



~2pp,+ - 2p') 



AfTO7v(pi+ - 4p) 



+TO^?7(2ry-5pi+)-2Af2+<^ 
+ £^[(p^-2p^ + 2p,.p), + 2A/2(2/ 
- 2pi+p) + 2T]m%{2pi+p'^ - 277/ + 

"^1+ J J 
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"1+ 

dadf3a lt^{<;;T) 



1 



W /^2 



2to2+ 



-ppi+-2p2)J _c;"(<,;T) (^^ 
p-2p^) + (?7TO7v - Af ) 



'■N 



2ml+ 



,2 

1+ 



(2p- 3pn 



~2ri)mN + 2M 



'1+ 



+TOjv[pi+?7 + + p^ipi+ - 2?7)]) 

2/p2^,^rn5^- 



(C16) 



= - 

vi2ai 



B"2 /.I 

-^1 + 



127r2 
2ppi+ - 2p2 



TO 



(?;r) (<;-2A// 
2 



1+ 
2 



~2AfTO7v(4p - pi+ - I877) + 277TOjv(5pi+ 
-6p - 57?) + ^[4TO^(p2+ - 2/ + 2pp^+) 



-1+ 



X (2A//2 - cr) + 36p?+,7AfTO^ + 477mj^ 
(77[p2^ -6pi+p-2p2]_4pi+p2^] 



32pf+phfm% 



m 



1+ 



(C17) 



127r2 



dadpa { Ci (^;T) 



'TV 



(2pi+P 



'1+ 



-2p2^ - 5/) 
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(M - mNv) 



lOM 



-2mN{3pi+ +r] + p) + ^^{M[2pi+p - bp^ 
-2pl+] + mw[7/(4pi+p + p2 + 4p2^) + 4pi+/]) 
\Qm%pl+p'^r|^ 1 



TO 



1 + 
2 ^2 



(C18) 



where (;B{M^,ml+,a,l3,pi+,-in%), p = I - rj, r] = 
a + af3pi+ . 

At T = 0, using the definitions in Table HI the nucleon 
Faddeev equation yields 



TOAr(T = 0) = 1.09 GeV, 



(C19) 



which is just 4% less than the value produced by the 
static approximation employed in Ref. [27[ . The Faddeev 
amplitude is 



s ai 
0.85 0.41 



02 

-0.34 



(C20) 



which may be compared with the static approximation 
resuh [13: (3,01,02) = (0.88,0.47,-0.078). Plainly once 
again, Eqs. (|C8p . (|C9p represent an improvement over 
previous definitions of the static approximation because 
they eliminate one parameter yet produce essentially the 



same mass. 



Here, too, we emphasise that Eqs. (jCSp . (|U9p describe 
the nucleon's "dressed-quark-core" . They should and do, 
therefore, produce a mass that lies above that quoted 
empirically. The inclusion of resonant contributions to 
the kernel leads typically to a 0.15 GeV reduction in 
the bound-state's mass [l^, H^]. Notably, the mass in 
Eq. (IC19[) is within 12% of the undressed nucleon's mass 
in Ref. 11J|. Moreover, TOjv/mp = 1.17, which compares 
favourably with the experimental value of 1.21. 



Appendix D: Euclidean Conventions 



Our T = Euclidean conventions are specified here. 



(Dl) 



{7m' li'} = 2 (5^1. ; 7^ = 7m ; ctm!- = 9 [7m' 7i^] ; (D2) 

tr [757M7i'7p7cr] = -4 e^^pg- , £1234 = 1 ■ (D3) 
A positive energy spinor satisfies 

u{P, s) (i7 • P + M) = = (i7 • P + M) u{P, s) , (D4) 
where s = ± is the spin label. It is normalised: 

u(P, s) u{P, s)=2M, (D5) 



and may be expressed explicitly: 



i(P, s) = VM - i'E g.p 



(D6) 



M - i 'E 



with £ = ?VP2 + M2, 




(D7) 



For the frce-particle spinor, ?2(P, s) — u{P, 5)^74. 

The spinor can be used to construct a positive energy 
projection operator: 



^+^^) ^= ^ E "(^' ^) = ^ ^-'^ -P + M). 

(D8) 



s=± 



A negative energy spinor satisfies 

v{P, s) {i-f-P- M) = = (n • P - M) v{P, s) , (D9) 

and possesses properties and satisfies constraints ob- 
tained via obvious analogy with m(P, s). 

A charge-conjugated Bethe-Salpeter amplitude is ob- 
tained via 



T{k;P) ^C'Ti^kiPf C , 



(DIO) 



where "T" denotes a transposing of all matrix indices and 
C = 7274 is the charge conjugation matrix, = — C. 
We note that 



(Dll) 



In describing the decuplet A-baryon we employ a 
Rarita-Schwinger spinor to represent a covariant spin- 
3/2 field. The positive energy spinor is defined by the 
following equations: 

{ij-P + M)Uf,{P;r)^0, (D12a) 
7^M^(P;r) = 0, (D12b) 
PpM^(P;r) = 0, (D12c) 

where r — —3/2,-1/2,1/2,3/2. It is normalised: 

Sp(P;r')M^(P;r) = 2M, (D13) 
and satisfies a completeness relation 

3/2 

2M 



a/2 

— J2 u^iP;r)uAP;r)^A+iP)R^,, (D14) 



-3/2 



where 



1 2 - - 1 - 

Rt,u = V-'^D - -^It^lu + -Pt,PuI\:, - T-^Pi^lv - Pvly\ , 

(D15) 

with P^ = —1, which is very useful in simplifying the 
Faddeev equation for a positive energy decuplet state. 
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